Evolution of rarefaction pulses into vortex rings. 
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The two-dimensional solitary waves of the Gross-Pitaevskii equation in the Kadomtsev- 
Petviashvili limit are unstable with respect to three-dimensional perturbations. We elucidate the 
stages in the evolution of such solutions subject to perturbations perpendicular to the direction of 
motion. Depending on the energy (momentum) and the wavelength of the perturbation different 
types of three-dimensional solutions emerge. In particular, we present new periodic solutions hav- 
ing very small energy and momentum per period. These solutions also become unstable and this 
secondary instability leads to vortex ring nucleation. 
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Considerable interest is attached to determining the 
entire solitary wave sequences of solutions of the Gross- 
Pitaevskii (GP) model 0] because they define possible 
states that can be excited in a Bose condensate. Jones 
and Roberts 12] determined the entire sequence of soli- 
tary solutions numerically for the GP model — 2i'0t — 
V'^ijj -I- (1 — |'!/'P)'0i where we use dimensionless variables 
such that the unit of length corresponds to the healing 
length a, the speed of sound is c = l/-\/2, and the density 
at infinity is poo = 1- The solitary wave solutions satisfy 
two conditions: (1) the disturbance associated with the 
wave vanishes at large distances: V' ~* Ij |x| — > oo and 
(2) they preserve their form as they propagate with a di- 
mensionless velocity U, so that ^{x, y, z, t) = iIj{x' , y, z), 
x' = X — Ut, in three dimensions (3D) and il)(x^y,t) = 
^(x', y), in two dimensions (2D), so that the solitary wave 
solutions satisfy 



2i(7^ - v^^ + (1 - iv-nv-. 



Jones and Roberts calculated the energy £ and mo- 
mentum V given by 

£ ^\ j \VM'' dV +\ j {l-\M''? dV (2) 

T' = ^JW-^)'^i'-i^-^)'^^*]dV, (3) 

and determined the location of the sequence on the V£- 
plane. 

In three dimensions they found two branches meeting 
at a cusp where V and £ assume their minimum values, 
VmXn and £min- hs V ~^ OO OTL cach branch, £ — > cx). On 
the lower branch the solutions are asymptotic to large 
vortex rings. 

As £ and V decrease from infinity along the lower 
branch, the solutions begin to lose their similarity to large 
vortex rings. Eventually, for a momentum Vq slightly 



greater than Vmim they lose their vorticity (-0 loses its 
zero) , and thereafter the solitary solutions may better be 
described as 'rarefaction waves'. The upper branch con- 
sists entirely of these and, as T' — » oo on this branch, the 
solutions asymptotically approach the rational soliton so- 
lution of the Kadomtsev-Pctviashvili Type I (KPI) equa- 
tion S and are unstable. In 2D the family of the soli- 
tary wave solutions are represented by two point vortices 
\i U < 0.4. As the velocity increases the wave loses its 
vorticity and becomes a rarefaction pulse. As f/ — > l/\/2 
both the energy, £ and momentum V per unit length ap- 
proach zero and the solutions asymptotically approach 
the 2D rational soliton solution of KPI. 

Jones and Roberts [|| derived the KPI equation us- 
ing an asymptotic expansion in the parameter e^ ~ 
2(1 — \/2[/), which is small when U approaches the speed 
of sound. They sought solutions of the form ifj = f + ig, 
where / = 1 + e^/i + e''/2 + • • •, 5 = e.9i 



e^.92 



(1) U ^ l/\/2 + e^C/i + • • •. The independent variables were 



stretched, so that ^ — ex' , rj = e^y, and ( = e'^z. By sub- 
stituting these expressions into (|l|) and considering real 
and imaginary parts at the leading and first orders in e, 
they determined that gi satisfies the KPI equation: 






Vf,^5i 






2 di^ x/2 V aC / . 



0, (4) 



and /i = dgi/V2d^ - gl, Ui = -l/2%/2. In 2D the 
equation (H) has a closed form solution B, so that the 
asymptotic solution of the GP equation in the KPI limit 

is 



i/> = 1-i- 



2V2a;' 



x" + e2y2 + 3/2e2 



g2y2 



3/2e2 



(5) 



which we have written in the original variables. 

It was shown by Kuznetsov and Turytsin B that the 
2D KPI soliton is stable to 2D but unstable to 3D per- 
turbations. The linear stability analysis of 2D solitary 



solution of the GP equation subject to long wavelength 
infinitesimal perturbations was done by Kuznetsov and 
Rasniussen g. They demonstrated that all long wave- 
length antisymmetric modes are stable and all long wave- 
length symmetric modes are unstable. In particular they 
showed that the growth rate of symmetric perturbations, 
a, is given by a^ = —£k'^/{dV/dU) > 0, as the wavenum- 
ber fc — > 0. The maximum growth rate of instability and 
the instability region of 2D solitary solutions were found 
in M by solving the linear stability problem. Through 
numerical integration of the GP equation it was shown 
that as perturbations grow to finite amplitude the vortex 
lines reconnect to produce a sequence of almost circular 
vortex rings. Senatorski and Infeld 1^] numerically in- 
tegrated the KPI equation to study the fate of 2D KPI 
solitons subject to 3D perturbations. They determined 
that 2D KPI solitons evolve into 3D KPI solitons which 
are also unstable. 

The goal of this Letter is to elucidate the fate of the 
2D rarefaction pulse in the KPI limit of the GP model 
subject to 3D perturbations. We discovered that such so- 
lutions may evolve into vortex rings and this establishes 
a new mechanism of vortex nucleation. We found that 
this mechanism can operate in different ways. The inter- 
mediate states may involve periodic solutions consisting 
of interacting 3D rarefaction pulses that belong to the 
lower branch of the Jones-Roberts cusp with V <Vo 

We have performed direct numerical simulations using 
the numerical method described in M . We solve the GP 
equation in the reference frame moving with the veloc- 
ity U[ chosen in such a way that the main disturbance is 
kept within the computational box: 



£ = V/V2 = 87re£/3. 



(8) 



-'l--<-'^ 



(i-ivnv- 



(6) 



In these computations we follow the evolution of the 
asymptotic solution (JSJ) extended along the z— axis and 
moving in the a;— direction. The dimensions of the com- 
putational box are Dx = 60, Dy — 60, Dz — 180. The 
ccy— faces of the box are open, to allow sound waves 
to escape; this is achieved numerically by applying the 
Raymond-Kuo technique ||l^. The z ~ and z = D^ 
sides are reflective. 

The soliton (||) was perturbed along the z— axis, so 
that at t = 



X —* X + 0.1 cos(fc2;). 



(7) 



We choose k so that N periods of this perturbation fit 
exactly into the D^ dimension of the box. There are two 
main parameters of the problem that determine the fi- 
nal outcome of the instability: e, which determines the 
configuration, energy, and momentum of the initial field, 
and the wavelength of the perturbation i — D^/N. It 
can be easily shown using (bh that the energy (0) and 
the momentum (H) of our initial field per wavelength of 
perturbation are given by 



First, we consider the evolution of the KPI solitary so- 
lution subject to large wavelength perturbations i — 
20, 30, 60 and e = 0.5. Fig. 1 illustrates the appearance of 
vortex rings through contour plots of the cross-section of 
the solution in the xz— plane with y — Dy/2 for £ ~ 20. 
The last panel shows the isosurface jf/'P = 0.2. 

FIG. 1. The contour plot of the density field of the 
cross-section of solutions of the GP equation. The time snap- 
shots show the evolution of the KPI solution (Bl) of the GP 
equation in the a;^— plane with y = Dy/2. In (Bl) we took 
e = 0.5 and the wavelength of the initial perturbation (R) is 
^ = 20. The last panel shows the isosurface \'ip\'^ — 0.2. The 
solutions starting with the third panel posess vorticity and 
evolve into equally spaced vortex rings. 
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According to the time snapshots of Fig. 1 the solution 
evolves directly into a set of vortex rings; other axisym- 
metric 3D solutions including 3D KPI solitons, are not 
involved. Moreover, exactly one vortex ring is generated 
for each wavelength of the perturbation. These vortex 
rings are distanced i healing lengths apart and have radii 
much smaller than £ (see Table 1). They therefore inter- 
act only weakly with each other. The energy and mo- 
mentum (g) of one period of the perturbation is used to 
create one vortex ring, the extra energy and momentum 
escapes and is carried away by sound waves (phonons). 
The vortex rings are aligned and propagate together with 
the same velocity. This arrangement of vortex rings is it- 
self unstable and cannot last forever. 

Similar calculations were done for i = 30 and £ = 60. 
The results are summarized in Table 1 which gives the en- 
ergy and momentum per wavelength of the perturbation 
of the initial field and the energy, momentum, velocity, 
and radius of the resulting vortex ring. 

Table 1. 



£ 


t^init 


Pinit 


C-ring 


Turing 


L/ring 


-fXring 


60 


251 


355 


99 


162 


0.45 


2.7 


30 


126 


178 


86 


132 


0.49 


2.35 


20 


84 


120 


71 


102 


0.53 


1.9 



Next we explore the evolution of the KPI limit soli- 
tary waves of the GP model subject to small wavelength 
perturbations. The effect of the decrease in the wave- 
length of the perturbation is twofold: 1) the energy and 
the momentum (@) per wavelength of the initial field 
is decreased, therefore leaving less energy available for 
creating a new entity and 2) these entities are in close 
proximity to each other so they strongly interact. Fig. 2 
plots the cross-sections of the solutions for e = 0.5 and 
£ = 60, 30, 15, 7.5, 3.75. The first two panels (Fig. 2a and 
2b) illustrate the vortex nucleation discussed earlier. To 
the best of our knowledge the periodic solutions shown on 
Fig. 2c-e are unknown in the literature on the GP model 
or the nonlinear Schrodinger equation. The interesting 
feature of these solutions is that they lack a vorticity and 
have small energy and momentum per period that tend 
to zero as ^ — > 0. These solutions can be understood as 
periodic pulse trains composed of the rarefaction pulses 
positioned on the lower branch of the Jones-Roberts se- 
quence with V < Vq. The interaction between adjacent 
pulses reduces the total energy per period. The analysis 
of these and other properties of periodic pulse trains com- 
posed from the solitary waves of nonintegrable evolution 
equations can be found in pl[. 



ilarly to the periodic solutions made of aligned vortex 
rings, the periodic rarefaction solutions become unsta- 
ble and we followed the development of this instabil- 
ity. Fig. 3a shows the contour plots of the density for 
a cross-section of the field which evolves from the peri- 
odic rarefaction solution with ^ = 15 to produce total of 
four rings. Therefore, to create each ring the energy and 
momentum of several periods of the rarefaction solution 
were used. The wavelength of the secondary instabil- 
ity that destroyed the periodic solution is approximately 
57. Similar calculations were done for the even shorter 
perturbation wavelength i = 7.5; see Fig. 3b. The wave- 
length of the secondary instability is approximately 29 
resulting in the appearance of six rings. The reason for 
the apparent nonuniformity of the nucleated rings is that 
Dz is not an exact multiple of these wavelengths. 

FIG. 3. The contour plot of the density field of the 
cross-section of solutions of the GP equation. The time 
snapshots show the different stages in the evolution of the 
KPI-limit solution (o) of the GP equation with e = 0.5 in 
the Tz— plane with y = Dy/2. The wavelengths of the initial 
perturbation ([^) are £ — 15 (a) and £ — 7.5 (b). 



FIG. 2. The contour plot of the density field of the 
cross-section of solutions of the GP equation. The time 
snapshots show the different stages in the evolution of the 
KPI-limit solution (0) of the GP equation in the a:z— plane 
with y — Dy/2. In (^ we took e = 0.5 and the wavelengths of 
the initial perturbation (Q) were ^ = 60 (a), ^ = 30 (b),£ = 15 
{c),£ — 7.5 (d), and £ — 3.75 (e). The contour plots are shown 
for |'0p at f = 149 for (a) and t = 127 for the rest of panels. 




Other findings are summarized in Table 2 which gives 
the values of the energy and momentum of the initial 
field and the resulting periodic rarefaction solution. 

Table 2. 



£ 


^init 


^init 




/ per 


15 


62 


90 


55 


78 


7.5 


31.5 


46.2 


30.7 


44 


3.75 


15.5 


23 


15 


21 



These periodic solutions can execute standing wave os- 
cillations of decreasing amplitude and period £/2. Sim- 





f=4i69 



t=554 



f = 4i69 



t = 554 



Finally, we consider the evolution of the small energy 
and momentum KPI limit solitary waves. In these com- 
putations we use e = 0.3 and £ = 60, 30, 20, 15. For 
£ = 60 the KPI solution follows the scenario of vortex nu- 
cleation and directly evolves into 3 vortex rings of small 
radii. For smaller wavelengths £ = 30, 20, 15 the KPI 
solution initially evolves into oscillating periodic rarefac- 
tion pulses of decreasing y— extent. The energy and mo- 
mentum per period are apparently insufficient to allow 
them to evolve into rings and the necessary energy can- 
not be reduced through interactions when the putative 
solutions are separated by such large distances. These 
solutions break down into sound waves that carry off all 
energy and momentum; see Fig. 4. 



FIG. 4. Isosurfaces of the density field of solutions of the 
GP equation. The time snapshots show the different stages in 
the evolution of the KPI-limit solution (H) of the GP equation 
with e = 0.3. The wavelengths of the initial perturbation (M) 
is ^ = 15. The minimum density increases with time and ap- 
proaches unity as the solution breaks down into sound waves. 
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In summary, we studied the instability of the 2D KPI 
hmit solitary wave solution in the GP equation. The evo- 
lution of several types of solutions are considered. We 
were able to identify three different regimes of transi- 
tion depending on the initial energy and momentum of 
the KPI solution and on the wavelength of the initial 
perturbation. For large wavelengths, the initial solution 
immediately evolves into a periodic solution consisting of 
small equally spaced vortex rings with a period equal to 
the period of the initial perturbation. For shorter wave- 
lengths the KPI solution first evolves into a periodic so- 
lution consisting of 3D interacting rarefaction pulses that 
later break up into vortex rings under the influence of a 
secondary instability of a different wavelength. Finally, 
if the energy of the KPI solution is small the solution 
can break into sound waves after forming an oscillating 
periodic rarefaction pulse. 

Fig. 5 summarizes all calculations performed and the 
relationships of the different regimes studied. The initial 
states considered are represented by dots on the V£— 
plane, where £ and V are defined per wavelength of the 
perturbation. We plot the cusp determined by Jones and 
Roberts [0| for the family of the vortex rings and rarefac- 
tion pulses. The arrows show the way the initial state 
evolve. 

FIG. 5. Summary of numerical integration of (H) starting 
with the initial condition (Bl) with e = 0.5 and e = 0.3. The 
cusp corresponds to the 3D solitary wave solutions. Dashing 
of the upper branch indicates that this branch is unstable. 
(Vo,£-o) marks the point where the vorticity disappears and 
point (Pminjfmin) gives the position of the lowest momen- 
tum-energy state of the 3D solitary solution (see discussion 
in the text). The line from the origin to ("Pmin, i^min) cor- 
responds to the family of periodic rarefaction pulses. Dots 
indicate the position of the initial states, the wavelength of 
the perturbation I is given next to each initial state; arrows 
show the evolution of these solutions, and the crosses cor- 
respond to the final state before the onset of the secondary 
instability. 
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